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INTRODUCTION 


Many  precise  satellite  attitude  estimation  systems  utilize  some  type  of  on-board  star  sensing 
device.  These  devices  provide  very  accurate  measurements  of  stellar  images  which  can  be  used  with 
the  inertial  coordinates  of  the  observed  star  to  obtain  vehicle  orientation  information.  A  key  step 
in  this  is  the  identification  of  the  imaged  stars. 

The  star  sensor  of  interest  here  can  be  thought  of  as  an  optical  device  which  consists  of  a  number 
of  light-sensitive  slits  on  an  otherwise  opaque  negative  image  plane.  As  the  vehicle  rotates,  images 
of  stars  cross  the  slits  and  the  precise  image  plane  coordinates  and  time  of  each  transit  are  recorded. 
Due  to  the  nature  of  these  sensors,  there  will  generally  be  no  simultaneous  hits  by  multiple  stars. 

This  nonsimultaneity  of  acquisition  compounds  the  star  identification  problem.  Traditionally,  in  this 
type  of  situation,  a  motion  model  for  the  vehicle  is  integrated  to  map  ail  of  the  star  hits  into  an 
estimated  celestial  coordinate  frame.  Then,  star  identification  is  attempted  by  either  a  direct  match 
(overlay  of  imaged  stars  and  catalog  stars)  or  by  comparisons  of  angular  separations  among  catalog 
stars  and  imaged  stars.  There  are  two  major  difficulties  with  this  mapping  technique.  First,  errors  in 
initial  orientation  can  cause  the  estimated  celestial  frame  to  be  substantially  displaced  from  the 
true  catalog  frame.  Second,  errors  in  the  rotational  motion  model  effectively  cause  distortion  in 
the  mapped  star  field.  Either  of  these  problems  can  cause  failure  of  a  direct  match  algorithm. 

Angular  separation  algorithms  are  sensitive  to  distortions  but  are  affected  by  initial  errors  only  in 
that  larger  subsets  of  the  reference  catalog  must  be  used  for  comparisons.  This  situation  can  quickly 
give  rise  to  a  prohibitive  number  of  computations.  Separately  or  in  combination,  these  errors  can 
easily  lead  to  nonidentification  of  the  star  pattern.  The  only  recourse  in  traditional  methods  is 
often  brute  force.  Parameters  of  the  motion  model  and/or  the  initial  orientation  are  systematically 
stepped  through  a  range  of  permissible  values  and  identification  is  repeatedly  attempted  using  one 
of  the  above  techniques. 

The  new  mapping  technique  described  here  provides  an  alternative  approach  where  star  hits 
are  collapsed  to  an  instantaneous  sensor-fixed  frame.  Inherently  contained  in  this  technique  is 
information  that  can  be  used  to  improve  parameters  of  the  motion  model  and  thus  decrease  distor¬ 
tion  effects.  A  two-stage  star  identification  process  then  allows  deterministic  computation  of  the 
absolute  vehicle  attitude  at  the  collapse  time.  This  in  turn,  decreases  computational  problems  due  to 
errors  in  initial  conditions. 


DEVELOPMENT  OF  EQUATIONS 


The  differential  equation  for  the  direction  cosine  matrix  that  orients  the  vehicle  fixed  refer¬ 
ence  frame  { g }  relative  to  inertial  space  {n }  is  given  by  Reference  1. 

(CJN(t)J  -  {u*t)]T  (GN(t)J  (1) 
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where 


0  «3(t)  -w2(0 

-w3(t)  0  w,(t) 
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w2(t)  ~Wj(t)  0 
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Wj(t),  (i  *  1 , 2, 3)  are  the  {g }  components  of  the  inertial  angular  velocity.  This  equation  is  purely 
kinematic.  Dynamics  enter  only  through  the  determination  of  the  angular  velocity  history.  It  is 
assumed  in  this  development  that  good,  dense,  but  discrete  samples  of  cJi(t)  are  available  from  an 
on-board  rate  gyro  system.  Under  this  assumption  the  entire  problem  becomes  kinematic  in  nature 
with  the  motion  characterized  by  Equation  (1),  which  contains  no  explicit  model  errors.  In  this 
case. 


«j(t)  B  WniW  “  Bj»  (*  ®  !>  2, 3) 

where  ojmi(t)  are  the  measured  {§}  components  of  the  angular  velocity  and  Bj  are  the  components 
of  the  rate  gyro  bias  vector. 

It  can  be  shown  that  (1)  has  a  state  transition  matrix  form  of  solution  given  by  Reference  2. 
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where 
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(3) 

with 

Wto.toH-m 

(4) 

At  a  given  star  hit  acquisition  time,  t,  the  inertial  and  sensor-fixed  direction  cosines  of  the  line  of 
sight  to  the  star  are  related  by 

K) 

-  [SG]  (GN(t)J  /  >  (5) 

(  S/ 

where  Sj|(t),  (i  *  1, 2, 3)  are  measured  or  “sensor-observed”  direction  cosines  to  the  star, 

figj,  (i  =  1, 2, 3)  are  the  unknown  but  time-independent  inertial  components  of  the  line  of  sight  to 
the  star  (the  “celestial  coordinates"),  and  [SG]  is  the  constant  rotation  matrix  which  orients  the 
sensor-fixed  reference  frame,  {£ } ,  relative  to  {g ) . 
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The  effective  sensor-fixed  direction  cosines  at  some  arbitrary  epoch  (say  to)  of  the  star  that  was 
actually  acquired  at  t  are 


f  yv) 

<  yv>= 

[  y*o>) 


[SGI  [GN(t0)J 


(6) 


Eliminating  the  unknown  inertial  direction  cosines  from  (6)  by  transposing  (S)  and  substituting 
into  (6)  gives 
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From  (2),  [GN(to)]  [GS(t)]T  can  be  identified  as  [«t>(t,  to)]T  so  that  (7)  becomes 
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with  [$(t,  to)]  determined  by  numerically  integrating  Equation  (3).  This  equation  provides  the 
basis  for  a  collapsing  technique  that  allows  all  star  hits  to  be  mapped  into  their  effective  sensor- 
fixed  positions  at  some  common  epoch  time  tp. 


Several  observations  can  be  made  concerning  star  identification  based  on  Equation  (8).  The 
integration  of  [4>(t,  tq)]  involves  only  the  angular  velocity  history  that  occurs  between  to  and  t. 
Furthermore,  since  the  differential  equation  for  [$>(t,  to))  is  kinematic,  there  are  no  inherent 
model  errors  in  Equation  (8).  The  only  errors  that  are  introduced  are  uncertainties  in  the  rate 
gyro  biases,  Bj,  and  any  residual  instrument  noise  that  may  be  present  in  the  measured  samples 
of  angular  velocity,  <om  r  In  the  absence  of  these  errors,  all  hits  made  by  a  given  star  as  it  crosses  the 
various  slits  in  the  image  plane  collapse  to  a  single  point  at  time  to  when  processed  by  Equation  (8). 
If  these  errors  are  present  but  reasonable,  multiple  hits  made  by  individual  stars  still  collapse  to 
distinct  clusters  of  points.  Consequently,  information  is  present  (namely,  the  dispersions  of  the 
clusters)  that  can  be  used  in  a  least  squares  adjustment  of  biases  without  a  priori  knowledge  of 
the  identities  of  the  imaged  stars.  With  unknown  biases  removed  in  this  fashion,  any  remaining 
dispersions  in  the  clusters  are  random  under  the  assumption  that  instrument  noise  in  the  gyros  is 
white.  The  best  estimates  of  the  true  image  plane  positions  of  the  stars,  if  they  had  been  ail  imaged 
at  tQ,  are  just  the  weighted  centroids  of  each  cluster.  The  net  result  of  this  is  an  instantaneous 
'‘photograph”  of  the  star  field  with  minimal  distortions  due  to  the  motion  model. 


Explicit  reference  to  vehicle  orientation  relative  to  inertial  space,  [GN(t)] ,  has  been 
eliminated  from  (8).  Clearly,  an  estimate  of  this  information  is  required  at  some  point  in  the  star 
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identification  process.  Since  the  result  of  the  collapsing  and  bias  adjustment  process  is  essentially 
a  distortion-free  instantaneous  photograph  of  the  star  pattern,  tentative  identification  of  a  pair 
of  stars  can  be  made  by  angular  separation  comparisons  among  imaged  and  catalog  stars.  The 
estimate  of  initial  orientation  is  now  a  problem  only  in  that  it  must  be  used  to  determine  the  proper 
subset  of  the  star  catalog  to  be  used  for  these  comparisons.  A  large  uncertainty  in  orientation  means 
a  large  subset  of  the  catalog  must  be  used.  However,  since  identification  of  only  a  single  pair  of 
stan  is  betog  undertaken  at  this  point,  the  computational  burden  is  minimized.  Once  two  stars  are 
tentatively  identified  by  this  technique,  a  deterministic  orientation  of  the  vehicle  at  to  can  be 
computed  as  follows.  Construct  both  sensor-fixed  and  inertial  components  of  a  right-handed  set  of 
orthonormal  unit  vectors  according  to 


“lsils  *lsIJI  +  *ls2*2  +  £ls3f3  =  £lni*l  +  £ln2*2  +  2ln3*3 


(9a) 


*2  *  £2*j  *1  +  £2»2  *2  +  C2s3  *3  =  C2nj  ®1  +  £2n2  ®2  +  22n3  ™3 
n2  *^1xS2/lK1Xfi2! 


(9b) 

(9c) 


u3sUjXt/2  (9d) 

where  and  $2  are  the  unit  vectors  to  stars  1  and  2,  respectively.  This  system  of  equations  can  be 
written  as 


<«}-  (A]{*}-  [B]{»}  (10) 

where  [A]  and  (B]  are  the{  s }  and  { h  }  components  of  the  unit  vectors. 

The  unknown  orientation  of  { s  }  relative  to  {  «  }  is  given  by 

U}  =  [SN(t0)J  { n  }  (11) 

This  expression  can  be  combined  with  (10)  and  solved  for  [SN]  to  yield 

[SN(t0)l  =  [A]t  [B]  (12) 

All  candidate  stars  can  now  be  mapped  into  the  tq  image  space  where  the  star  pattern  is 
accepted  or  rejected  based  on  a  direct  match  scheme.  If  the  pattern  is  not  accepted  as  valid,  two 
new  stars  are  tentatively  identified  by  angular  separation  comparisons  and  the  process  is  repeated 
from  that  point. 


CONCLUSIONS 


The  development  outlined  herein  is  based  on  the  availability  of  observations  of  the  vehicle 
angular  velocity  vector.  This  general  star  identification  process,  however,  is  not  restricted  to  that 
situation.  If  angular  velocity  measurements  are  not  available,  a  dynamic  model  for  the  vehicle  mo¬ 
tion  must  be  incorporated.  The  integration  of  Equation  (3)  then  depends  on  the  resulting  dynamic 
angular  velocity  history.  The  dispersions  within  clusters  of  star  hits  after  collapsing  still  contain 


information  that  can  be  used  to  improve  parameters  of  the  dynamic  model.  Implementation  of  the 
remainder  of  the  process  remains  essentially  unchanged. 

This  technique  where  traditional  approaches  are  combined  and  modified  by  the  generation  of 
an  instantaneous  photograph  is  a  viable  approach  to  star  identification.  In  cases  where  the  collapsing 
is  rigorous  due  to  the  availability  of  good  angular  velocity  measurements  and  where  instrument 
biases  can  be  recovered  because  of  multiple  hits  by  single  stars,  this  approach  is  a  particularly 
attractive  alternative. 
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